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Abstract 

We study numerically the so-called fundamental modular region A, a region free of 



I Gribov copies, in the minimal Landau gauge for pure SU{2) lattice gauge theory. To this 

CnI ■ end we evaluate the influence of Gribov copies on several quantities — such as the smallest 

eigenvalue of the Faddeev-Popov matrix, the third and the fourth derivatives of the min- 
imizing function, and the so-called horizon function — which are used to characterize the 
I region A. Simulations are done at four different values of the coupling: (3 = 0, 0.8, 1.6, 2.7 , 

and for volumes up to 16^. We find that typical (thermalized and gauge- fixed) configura- 
, tions, including those belonging to the region A, lie very close to the Gribov horizon d^l, 

and are characterized, in the limit of large lattice volume, by a negative-definite horizon 
tensor. 



.J 1 Introduction 

a ^ Gauge theories, being invariant under local gauge transformations, are systems with redundant 
dynamical variables, which do not represent true dynamical degrees of freedom. The objects 
of interest are not the gauge fields themselves, but rather the classes (orbits) of gauge-related 
fields. The elimination of such redundant gauge degrees of freedom is essential for understanding 
and extracting physical information from these theories. This is usually done by introducing a 
gauge-fixing condition which determines a representative gauge field on each orbit. In reference 
|1[| Gribov showed that Coulomb and Landau gauge-fixing conditions do not fix the gauge fields 
uniquely, i.e. there are many gauge-equivalent configurations satisfying the Coulomb or Landau 
transversality condition. These Gribov copies do not affect perturbative calculations, but their 
elimination could play a crucial role for non-perturbative features of gauge theories. Gribov's 
result was generalized by Singer M to the case of a generic compact semi-simple non-abelian Lie 
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group. In this work the author considered continuous gauge fixing, and boundary conditions 
(i.e. value of the gauge fields at infinity) such that the euclidean space-time is compactified to 
the four-dimensional sphere 5^. A similar analysis was done by Killingback for the case 
of periodic boundary conditions (i.e. the four-dimensional torus T^), which are the boundary 
conditions usually employed in lattice gauge theory. 

A possible solution to the problem of Gribov copies is to restrict the functional integral to a 
subset of the gauge-field space, the fundamental modular region A, which is the set of absolute 
minima of a Morse function on the gauge orbits ^. In the so-called minimal Landau gauge 
this Morse function is defined as 

SA[g] = M(^Mr f d'x{[4^^Y{x)y . (1) 

It has been proven §], 0) @] that, with this choice of the Morse function, every orbit intersects 
the interior of the fundamental modular region once and only once, i.e. in the interior of A the 
absolute minima are non-degenerate and there are no Gribov copies. However, on the boundary 
of the fundamental modular region there are degenerate absolute minima, and only when they 
have been identified can we obtain a region truly free of Gribov copies [^, 

The problem of Gribov copies is also present in the lattice regularization of gauge theories 
IP, |10|- Although this formulation does not require gauge fixing, due to asymptotic freedom, 



the continuum limit is the weak-coupling limit, and a weak-coupling expansion requires gauge 
fixing. Thus, gauge- dependent quantities are usually introduced on the lattice, and Gribov 
copies have to be taken into account in lattice gauge theory as well. 

A fundamental modular region can be defined also on the lattice. In this case, for the 
minimal Landau gauge, we consider the absolute minima of the functional^ 

^uig] = ^EETr - g{x) U,{x) g^x + e,)] , (2) 

which is the lattice analogue of the Morse function SA[g] used in the continuum [see eq. (|l|)]. 
Since the gauge orbit is a compact manifold on a finite lattice, this functional is bounded 
and the existence of an absolute minimum for Su follows immediately. Let us notice that the 
functional Su[g] can be rewritten |T2[ as the quadratic formQ 



^u[9] = ^T.9ix) ■{Vug){x), (3) 

where 

{Vug){x) = E [^9\x) - U^{x)g\x + e^) - Ul{x - ae^) g\x - e^)] (4) 

is the gauge-covariant Laplacian. The minimization of a quadratic form is a standard and 
simple problem if the variables are elements of a linear space In our case, however, these 



^ This definition applies to SU{2) lattice gauge theory in 4 dimensions. We consider a standard Wilson 
action with periodic boundary conditions, and lattice volume V . For notations we refer to [ pj| . 
^ In equations (0) and (0) we use g to indicate four-dimensional unit vectors. 
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variables are SU{2) matrices (i.e. 4- dimensional unit vectors), and due to this non-linearity the 
numerical search for the absolute minimum becomes highly non-trivial. 

Most of the properties proved in the continuum for the fundamental modular region can be 
extended to the lattice case [|14|, ITSf. In particular, an explicit example of degenerate absolute 
minima on the boundary of A is given in reference . 

In this work we want to characterize the fundamental modular region by evaluating numer- 
ically diagnostic quantities (see Section ^ at relative and absolute minima. To date, relatively 
few studies |]TD|, [IT], |T^ have been done in order to analyze the influence of Gribov copies ( Gri- 
bov noise) on some lattice quantities. However, these numerical studies have never tried to 
characterize the "geometry" of the orbit space. We know general properties of the fundamental 
modular region, and we know particular cases in which this region can be studied analytically 
& i) 0) H) 0) HH 1131 • contrary, we do not know what happens in numerical simulations. 



This work aims at filling this gap, and at providing information about the "localization" in the 
gauge-field space of a typical thermalized gauge-fixed configuration. Preliminary results have 
been reported in [113]. 



2 Lattice Landau Gauge-Fixing Condition 

In this section we analyze in more detail the lattice Landau gauge-fixing condition. Let us 
recall that this gauge condition is imposed by minimizing the functional Su, defined in eq. (^, 
with respect to the variables {g{x)}, keeping the thermalized configuration {Ufj,{x)} fixed. 
We consider a one-parameter subgroup g^r; x) of SU{2) defined by 

g{T]x) =exp[ir7(2;) - ct] ; (5) 

here the parameter r is a real number, 7(x) is a three-dimensional real vector, and the com- 
ponents of a are the three Pauli matrices. Then the functional Ejj, defined in eq. (|^), can be 
regarded as a function of the parameter r, and its first derivative — with respect to r and at 
r = — is given by the well-known expression 

^'(0) = T^EEr(^) (V-A(^))"(x), (6) 

X a 

where 

{V-Ar{x)^Y.K(^)-Al{x-e,) (7) 

is the lattice divergence of the gluon field A^(a;) . If {t/^ (x)} is a stationary point of £{t) at 
r = then we have i^'(O) = for every {7"(x)}. This implies 

(V ■ Af {x) = (8) 

for any x and a, which is the lattice formulation of the usual Landau gauge-fixing condition in 
the continuum. 
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Still considering the one-parameter subgroup defined in eq. (^), we can evaluate the second 
derivative of S{t). Following references [|14[ |l^ one can check that 



1 

4V 



(9) 



^iV a,b 



6, c 



fi, a 
abc 



(10) 



where Aii/'{x,y) is the lattice Faddeev-Popov matrix and e"*'^ is the complete anti-symmetric 
tensor. It is clear that this second derivative is null if the vector 7"(a;) is constant, i.e. the 
Faddeev-Popov matrix has a trivial null eigenvalue corresponding to a constant eigenvector. 
From eq. (P) we obtain that, if {?7^ (x)} is a local minimum of S{t) at r = 0, then the matrix 
M-u is positive definite (in the subspace orthogonal to the space of constant vectors). This 
implies that any local minimum of the functional Su belongs to the region (see Figure jT]) 

n = {A : V ■ A = 0, Mu > 0} . (11) 

This region was introduced by Gribov ^ in the attempt of getting rid of spurious gauge copies. 
It is delimited by the so-called first Gribov horizon dfl, i.e. the set of configurations for which the 
smallest nontrivial eigenvalue Xmm of the Faddeev-Popov matrix is zero. Clearly the Gribov 
region Q includes the fundamental modular region A. However, there are points — the so- 
called singular boundary points P, ^ |l^, |15|, |17| — on the boundary of A which are also on the 
boundary of Q (see Figure p. A typical example are degenerate absolute minima of the Morse 
function for which the degeneracy is continuous. In fact, as said in the Introduction, all the 
degenerate absolute minima are found on the boundary of the fundamental modular region. In 
the case of a continuous degeneracy, the Faddeev-Popov operator must have a zero eigenvalue, 
i.e. these minima are also on the boundary of Q. An explicit example of singular boundary 
points on the lattice has been found by Zwanziger (see |jl^. Appendix E]). 

Finally, we can evaluate the third and fourth derivatives of S{t) with respect to r, at r = 0. 
Following again |14, |15| we obtain]^ 



£ (0) 



AV 



/i, a 



{i + e,) + ','(x) Al(x) 



V ^ 



7 [.X, 



r[x] 



(V ■ Af 



(12) 



and 



^""(0) 



l|FEE|[7V + e.)-f(x)]V» 

1 ^ I ^2, 



^Ejrix 

a;, J/ a,b 



(13) 



^ Our notation is slightly different from the notation used in references [Q and |T^. This explains the 
difference between the coefficients in formulae (|l2|) and ( p^ and the coefficients in the corresponding equations 
in references |l^ . 
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3 Characterization of the Fundamental Modular Region 



In order to characterize the "locahzation" of the gauge-fixed configuration in the hyperplane 
r of transverse configurations {V ■ A = 0; see Figure |I]), we consider the smallest nonzero 
eigenvalue Xmm of the Faddeev- Popov matrix M.u, and its corresponding eigenvector u;^j„(x), 
i.e. we solve the eigenvalue problem 



(14) 



As already said, the matrix has a trivial null eigenvalue corresponding to a constant 
eigenvector. Therefore, this equation has to be solved in the subspace orthogonal to constant 
vectors, namely the eigenvector uj^^^{x) should satisfy the relation J2x ^mini^) = 0- for 
the eigenvalue Xmin, we know that it is positive, since any local minimum of £u belongs to fl. 
Moreover, for the vacuum configuration U^{x) = IL (which also belongs to Q and to A), the 
Faddeev-Popov matrix is simply (minus) the lattice LaplacianQ. Therefore, in this case, the 
value of Xmin is given by 



XLap = 4sin2(7r/Ar) 



(15) 



where N is the lattice size. Finally, this eigenvalue goes to zero as the first Gribov horizon 
dfl is approached. So, the value of Xmin can be interpreted as a sort of distance between the 
minimum and dQ, and we may test whether its value is largest (in average) for the absolute 
minimum. 

After we have evaluated the eigenvector ci;^j„(x), we can set 7"(x) = ci;^j„(a;) in the one- 
parameter subgroup defined in eq. (|^). Then, if { U^{x) } is the configuration which minimizes 
the functional Su, we can study the behavior of the minimizing functional near the minimum 
using the gauge transformation generated by this one-parameter subgroup. More exactly, we 
can analyze its behavior along the "direction" of ti;^j„(x), i.e. along the direction for which the 
rate of increase of the functional £u is smallest. To this end we expand £{t) in powers of r 
around the minimum, i.e. around r = 0. We then write 



S(r) 



£{0) + 



E (0) + -E (0) + -E (0) 



+ . . . 



(16) 



where the derivatives of £^(r) [see equations (0), (0) dH)] evaluated with ^°-{x 
'^mm(^)' a^'i the eigenvector oj^^^Jyx) is normalized to one. We can now define the ratio 



^"'(0) 



" r(o)r"(o) ' 

which is independent of the scale r, and rewrite equation ( [T6| ) as 



:i7) 



E{t\ 



m + 



r2r(0) 



1 tE"\^) 

3 r(o) 



12r 



rr(0) 

r(o) 



+ 



This can be seen from equation (10) with U^Ax) = 1 and A^(a;) = 
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If we make the change of variables x = t £"' (0) / S" (0) , then it is clear that the shape of the 
minimizing function around its minimum is fixed by the value of the ratio r. As an example, 
in Figure we show the behavior of S{t) — S{0) for six different values of the ratio r. Let us 
notice that, for r > 8/3, there are Gribov copies of the minimum at r = 0; in particular, there 
is a maximum, which does not belong to the Gribov region Q, and a second minimum, which 
is an example of a Gribov copy inside the first Gribov horizon dQ. Thus, we expect a value of 
r smaller (in average) for the absolute minimum than for a generic relative minimum. Let us 
also notice that plots in Figure ^ are related to the bifurcation process described in reference 
In that case the author was following a stationary point of the minimizing function, while 
moving from inside to outside the region Q. Here, on the contrary, we sit at the minimum at 
r = and look at its surroundings for different values of the ratio r. We note that changing 
the value of r is equivalent to moving this minimum inside the region Q. 

If we set ■y°'{x) = u!'^^J^{x), and uj^^^{x) is normalized to 1, then from formulae and ([l^) 
we obtain A^m = 4:V S (0). In our simulations we evaluate the eigenvalue Xmin, the third and 
fourth derivatives of the minimizing function^, and the ratio r. 

Recently Zwanziger proposed |T^, ^ a modification of the SU (A^)-Yang-Mills action which 
effectively constrains the functional integral to the fundamental modular region A in minimal 
Landau gauge. This effective action is given by 

S,ff[U] = (3Sw[U]+aH[U] , (19) 

where (3 Sw[U] is the standard Wilson action, a is a new thermodynamic parameter and H[U] 
is the so-called horizon function defined as 

Hp] ^ ^j:H,a,,a[U]. (20) 

Here H^a^ubiU] is the horizon tensor given by 



and 



Ul{x) - Ul{x-ae,) 



+ e 



A'ix) + A'ix-ae, 



(22) 



With this new action the standard Yang-Mills theory is recovered (in the infinite- volume limit) 
only when the thermodynamic parameter a has the critical value fixed by the so-called 
horizon conditior^ 

-a-' = hm ^{H[U]). (23) 

For the vacuum configuration U^{x) = IL we obtain B'^^^[x) = 0, and therefore the horizon 
tensor is diagonal and equal to — 6^^6ab V. On the contrary, at the first Gribov horizon dQ, the 
term containing the inverse of the Faddeev-Popov matrix blows up, i.e. both the horizon tensor 



^ Also these derivatives are evaluated at r = and for 7° (a;) = a;J^jj„(a;) [see eq. ( p^ ) and (p^)], and 
multiplied by 4 V. 

^ Here (.) indicates an expectation value calculated in an ensemble Z{a) depending on the thermodynamic 
parameter a. For a definition of Z{a) see eq. (1.14)]. 
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and the horizon function are positive and become larger and larger as we approach d^l. In 
particular, points on the boundary of Q where the horizon function is infinite can be explicitly 
exhibited [|19| . However, the horizon function is not necessarily infinite for all configurations on 
the boundary of Q. In fact we can rewrite this function as fTl 



Hp] = r2^^Y 



E 



\ E E K^-) 

X fj, 



(24) 



where A„ and uini^) ^^^5 respectively, the eigenvalues and the corresponding normalized eigen 
vectors of the Faddeev- Popov matrix Aiu- Then, for the so-called degenerate gauge orbits flE 



one obtains that the absolute minimum is degenerate, namely it is on the common boundaries 
of A and Q, and that the term of the horizon function which is proportional to 1/Amm is of 
the indeterminate[] form 0/0. In the infinite- volume limit the horizon tensor is expected |T^ to 
be negative-definite inside the fundamental modular region A, and to vanish on its boundary. 
In the same limit, the measure should get concentrated [^] on a region where the horizon 
tensor per unit volume h^a,ub{U) is equal to zero. It is interesting to recall that the condition 
h{U) = 0, where h{U) is the horizon function per unit volume, is satisfied for all transverse 
configurations on a finite lattice with free boundary conditions . However, this result is not 
sufficient to make lattices of this kind free of Gribov copies [^. In order to test these conjec- 
tures, we evaluate the 12 eigenvalues of the horizon tensor (per unit volume) and the horizon 
function (also per unit volume). We also consider the average projection (per unit volume) of 



the B;j 



x) vectors on the eigenvector cj^ 



I.e. 



V 



12 V 



E 



E 



(x) uj" ■ 



x) 



(25) 



and the contribution V/Xmin to the horizon function h from this eigenvector [see formulae 
(p^) and (pSl) ]. Finally, we evaluate the largest and the smallest (both in absolute value) non- 

and the average over its non-diagonal 



diagonal elements of the matrix h^a,ub 
elements defined as0 

1 
132 



: 12 h 

E 



fiOjU b: 



1/2 



(26) 



4 Numerical Simulations 



In order to evaluate the horizon tensor, and the horizon function, we have to invert the Faddeev- 



Popov matrix A4u [see equations ( poD and (|2l|) ]. This matrix is rather large and sparse, and 
has an eigenvalue zero corresponding to constant modes. This means that it cannot be inverted 
directly. Nevertheless, this inversion can be done by using a standard conjugate- gradient (CG) 
method, provided that we work in the subspace orthogonal to constant vectors (see refences 



^ This result follows from the relatfon Bf^^{x) 



- [Dj^[U] )^, where -D^[C/] is the lattice gauge-covariant 
derivative, and the fact that degenerate gauge orbits have a nonzero solution for the equation D^[[/]a; — 0. 
For details see reference p5[ |. 

^ Here the index i (respectively j) stands for the pair of indices /i and a (respectively v and b). 
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TT| , |T8|, In particular, we have to impose that the source and the initial guess for the 

CG-method have zero constant mode. In our case the source is given by one of the (twelve) 
vectors B'^^{x) defined in eq. (^). Thus, before inverting the matrix A4u we have to impose 
the condition 

E B',ai^) = ■ (27) 

X 

The evaluation of the smallest nonzero eigenvalue of A1^^(x, y) can be done using the routine 
that inverts this matrix. To this end let us consider the vector 



) = EE(^^r)"'(^,2/)^'(2/), (28) 



m, we have 



where w'^{y) is a randomly chosen vector such that J2y ^'^{u) = 0. Then, in the limit of large 



= cA^™,a;l„(a;) + ... , (29) 
where c is an unknown constant. This implies 

for any a and x, and c<j^in(^) given by the vector uj'^{x) normalized to 1. There are of course 
more sophisticated algorithms to evaluate eigenvalues and eigenvectors of a matrix. However, 
for this case, this simple method converges very fast^ and is easy to be implemented, since it 
uses the same conjugate-gradient routine used to evaluate the horizon tensor. 
We consider, for each quantity, two different averages: 

• Average (indicated as "am") only on the supposed absolute minima.|^ 

• Average (indicated as "fc") only on the first gauge- fixed gauge copy generated for each 
thermalized configuration. This is the result that we would obtain if Gribov copies were 
not considered. 

If the result obtained from these two averages are systematically different for a given quantity, 
we can say that the existence of Gribov copies introduces a bias (Gribov noise) on the numerical 
evaluation of that quantity. 

The parameters used for the simulations can be found in |Tl|, Table 1]. However, in this 



paper, configurations at (3 = 1.6 and with lattice volume V = 24^ were not analyzed. Our runs 
were started with a randomly chosen configuration. Details about thermalization (using hybrid 
overrelaxation) and gauge-fixing (using stochastic overrelaxation) can be found in references 
|Tl|, Computations were performed on several IBM RS-6000/250-340 workstations at 

New York University. 



^ This tells us that the second smallest nonzero eigenvalue of J^i/'{x,y) is usually much larger than A 



mm • 



10 We refer to jl^, |l^ for details about the numerical search for a candidate for the absolute minimum of the 
minimizing function. 
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5 Results 



In Table |l| we report the data for the ratio Xmin/^Lap, the third and the fourth derivatives of the 
minimizing function, and the ratio r defined in eq. (|T^). For the third derivative we consider 
the absolute value since this quantity is negative for about 50% of the configurations. The 
fourth derivative, on the contrary, is always positive, with very few exceptions at = 2.7 and 
small lattice volumes. 

Our results are consistent with the conjectures discussed in Section ||: the value of Xmin 
(respectively r) is larger (smaller), in average, for the absolute minima (average "am"). On the 
contrary, Gribov noise is not observable for the third and fourth derivatives of the minimizing 
function. Thus the Gribov noise for the ratio r is entirely due to the Gribov noise of the 
eigenvalue Xmin- It is also interesting to observe that, for (3 in the strong-coupling regime, Xmin 
is down by a factor approximately 40 — 60 with respect to the lowest nonzero eigenvalue XLap 
of the negative of the lattice Laplacian. (At /3 = 2.7 this ratio is still about 6.5.) This indicates 
that typical (thermalized and gauge-fixed) configurations, including those that are absolute 
minima, lie very close to the Gribov horizon dQ, where Xmin = 0. Recall that the fundamental 
modular region A is included in the Gribov region Q. Consequently we conclude that a typical 
configuration in A lies close to its boundary dA, and where this boundary in turn lies close to 

dn. 

As for the quantity r, from Table |I] we observe that (where the statistics are good) r is 
quite small|3, even though typical configurations lie near the boundary of the Gribov region 
dfl, as indicated by the remarkable smallness of Xmin = 4:V S (0). At a generic point of dQ 
we have S (0) 7^ 0, so r = 00. However, at those points of dQ that are also points of dA, 
we have S (0) = (see references f^, |l4[ |l5|). The smallness of r for relative (as well as 



absolute) minima is consistent with the hypothesis that also for typical configurations that are 
relative minima the boundary dA lies close by. Note that whereas the gauge orbit is tangent 
to r (9 ■ A = 0) on a generic point of dQ, it is also tangent to dQ at the common points of dA 
and dn B |151. 

In Table ^ we report the data for the horizon function h, the contribution V/Xmin to the 
horizon function, and the minimizing function Su. Finally, in Table ^ we show the values for 
the smallest and the largest eigenvalues of the horizon tensor ha^^bv, and the average over its 
non-diagonal elements, as defined in eq. (p6|). 

From our data there is evidence of Gribov noise for h and for the eigenvalues of the horizon 
tensor: these quantities, in fact, are smaller at the absolute minima. Moreover, the horizon 
tensor is closer to a diagonal matrix at the absolute minima^. On the contrary there is no 
evidence of Gribov noise for the quantity V defined in eq. (P^. This suggests that the Gribov 
noise for V jXmin-, a quantity which contributes to the value of the horizon function, is due to 
the eigenvalue Xmin- Note that Ey — \ contributes to the horizon function h also [see formulae 
(0) and (p4D], and that Eu is (by definition) smallest at the absolute minimum. However, the 
Gribov noise for h is much larger than the Gribov noise of Eu and V jXmin-, i-e. it is probably 
due to Gribov noise for all the eigenvalues of the Faddeev-Popov matrix. 

As for the volume dependence, it seems that the horizon tensor becomes closer and closer 

We have found only in very few cases a value of r larger than 8/3, which indicates the existence of other 
stationary points near the minimum at r = (see Figure ^ . 

See data for the average of the non-diagonal elements in Table ^. 
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to a diagonal matrix as the lattice size increases (at fixed (3). On the contrary, the value of 
the horizon function seems almost independent of the volume, while V and V/\min decrease 
with V . It is also evident that, at large enough volume, the horizon tensor is negative-definite. 
However, for small lattice volumes we find positive eigenvalues of the horizon tensor. We have 
also found, again at small values of V , a few configurations with a positive value for the horizon 
function /i, even though h is always negative in average, and its value approaches —1 (the value 
on the vacuum) as (3 increases. 



6 Conclusions 

Our data show Gribov noise for quantities related to the Faddeev-Popov matrix, such as the 
eigenvalue Xmin and the horizon tensor. This is in agreement with the result obtained in 



reference |Tl]], in which Gribov noise has been observed for the ghost propagator. In all cases, 
the effect is small but clearly detectable for the values of (3 in the strong-coupling region. The 
fact that this noise is not observable at /? = 2.7 seems to us to be related only to the small 
volumes considered here. Of course this hypothesis should be checked numerically. This is, at 
the moment, beyond the limits of our computational resources. 

As for the "localization" in the gauge-field space of a typical (thermalized and gauge-fixed) 
configuration, the smallness of Xmin and r suggest that these configurations are always close to 
the common boundary of the Gribov region VL and the fundamental modular region A. This is 
the case for relative as well as for absolute minima of the minimizing function. Moreover, in 
the limit of large lattice volume, only the region 

Q = [A: AeQ,H^a,AU] < O} (31) 
seems to contribute to the evaluated expectation values. Thus our results support the conjec- 



tures |T5[ that, in the limit of large lattice volume, the horizon tensor is negative-definite inside 
the fundamental modular region, and that the measure should get concentrated on the common 
boundary of f2 and A. 
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Figure 1: The plane of the figure represents the hyperplane T of transverse configurations 
(V ■ A = 0). It contains the Gribov region Q defined in eq. (0), and the fundamental modular 
region A. Note that parts of the boundaries dQ and dA are in common. 
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Figure 2: Plot of S{t) — £{0) as a function of r for six different values of the ratio r [see 
equations (|T^ and (|T8|)]. For each value of r, the values for £ (0) and £ (0), and the scale on 
the y axis have been chosen "ad hoc" in order to emphasize the differences between the various 
cases. Recall that we are sitting at the minimum at r = 0. 
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Table 1: The ratio between the smallest nonzero eigenvalue Xmin of the Faddeev-Popov matrix 
and the smallest nonzero eigenvalue X^ap of (minus) the lattice Laplacian [see eq. (|15D], the 
absolute value of the third derivative and the fourth derivative of the minimizing function 
£u [see formulae (0) and ([T3|)], and the ratio r defined in eq. (|1^. The derivatives of the 
minimizing function are multiplied by 4 V, where V is the lattice volume. Two different types 
of statistics are considered: am = absolute minimum, and fc = first copy. Error bars (in 
brackets) are one standard deviation. 
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Table 2: The horizon function (per unit volume) h = H/V [see eq. (0)], the contribution 
V /Xmin to the horizon function h from the eigenvector uj'^^^{x) [see formulae (|^) and (|^)], and 
the minimizing function £u defined in eq. (0). Two different types of statistics are considered: 
am = absolute minimum, and fc = first copy. Error bars (in brackets) are one standard 
deviation. 
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Table 3: The smallest and the largest eigenvalues of the matrix hfj^a,ub = ^'^ H^a,ub/V [see 
eq. (piD], and the average of its non-diagonal elements [see eq. (^^]- Two different types of 
statistics are considered: am = absolute minimum, and fc = first copy. Error bars (in brackets) 
are one standard deviation. 
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